Magnetic Fields in Star-Forming Molecular Clouds III. 
Submillimeter Polarimetry of Intermediate Mass Cores and 

Filaments in Orion B 



Brenda C. Matthews 

, Department of Physics and Astronomy, McMaster University, Hamilton, ON, Canada L8S 4^1 

I matthewsOphysics .mcmaster . ca 

(N 

^ ■ Jason D. Fiege 

Canadian Institute for Theoretical Astrophysics, University of Toronto, Toronto, ON, Canada M5S 3H3 

' fiege@cita.utoronto.ca 

and 
T— I ■ 

> 

^ ■ Gerald Moriarty-Schieven 

^■f^ , National Research Council of Canada, Joint Astronomy Centre, Hilo, HI U.S.A. 96720 

g.moriarty-schievenSjach.hawaii . edu 

o ■ 

O 

Ph; abstract 

O ' Using the imaging polarimeter for the SubmiUimeter Common User Bolometric Array at the 

^ , James Clerk Maxwell Telescope, we have detected polarized thermal emission at 850 /im from 

^ ' dust toward three star-forming core systems in the Orion B molecular cloud: NGC 2071, NGC 

2024 and LBS 23N (HH 24). The polarization patterns are not indicative of those expected 
for magnetic fields dominated by a single field direction, and all exhibit diminished polarization 
, percentages toward the highest intensity peaks. NGC 2024 has the most organized polarization 

' pattern which is centered consistently along the length of a chain of 7 far-infrared sources. We 

^ have modeled NGC 2024 using a helical field geometry threading a curved filament and also as a 

magnetic field swept up by the ionization front of the expanding H II region. In the latter case, 
the field is bent by the dense ridge, which accounts for both the polarization pattern and existing 
measurements of the line-of-sight field strength toward the northern cores FIR 1 to FIR 4. The 
direction of the net magnetic field direction within NGC 2071 is perpendicular to the dominant 
outflow in that region. Despite evidence that line contamination exists in the 850 /im continuum, 
the levels of polarization measured indicate that the polarized emission is dominated by dust. 

Subject headings: ISM:clouds - ISM:magnetic fields - polarization - stars:fomiation - submillimeter 

1. Introduction the nearest giant molecular clouds and is an active 

site of low- to high-mass star formation. It was 
The Orion B (L1630) molecular cloud, at a dis- the first clouds to be systematically studied 

tancc of 415 pc (Anthony-Twarog 1982), is one of ^^^^^ ^^^^^ by Lada et al. (1991a), who found 



^Current address: Dept. of Astronomy, UC, Berkeley n i , i i . i . ;i i 

(bmatthews@astro.berkeleyedu) largest clumps, which together make up more 



that massive star formation takes place only in the 
five largest clumps, which together make up more 
than 50% of the mass of dense gas. We have cho- 
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sen three of these five clumps, NGC 2071IR (LBS 
8), NGC 2024 (LBS 33), and HH24 (LBS 23), for 
the current study. A fourth region, NGC 2068, 
contains a string of substantially smaller, fainter 
cores connected by weak dusty filaments (Mitchell 
et al. 2001). Polarimetry of this region will be 
presented in a forthcoming paper (Matthews & 
Wilson 2001, Paper IV in this series). 

Although the three regions have comparable 
gas masses, ranging from 230-460 M0 (Lada et 
al. 1991a), they have very different star forma- 
tion properties. NGC 2071IR lies four arcminutes 
north of the reflection nebula NGC 2071. This ex- 
tended submillimeter source consists of a cluster 
of at least nine embedded infrared stars (Walther 
et al. 1993) with a combined infrared luminosity 
of 520 Lq (Butner et al. 1990). The source IRS 
3 is thought to be the driving source of a mas- 
sive bipolar molecular outflow (Bally 1982; Snell 
et al. 1984; Moriarty-Schieven, Snell & Hughes 
1989; Eisloffel 2000). Houde et al. (2001) infer 
alignment between the outflow and its magnetic 
field by comparison of spectral lines of neutral and 
ionic species. Shocked molecular hydrogen (Bally 
& Lane 1982) and H2O mascrs (Genzel & Downes 
1979) are also seen towards this region, which is in 
a later evolutionary stage than NGC 2024 (Laun- 
hardt et al. 1996). Eisloffel (2000) has also docu- 
mented several other outflows in the region. Sub- 
millimeter maps of this region have been published 
by several authors (Mitchell et al. 2001; Motte 
et al. 2001; Johnstone et al. 2001). Based on a 
comparison between submillimeter continuum and 
CO(3-2) and HCO+ Hue data, Motte et al. (2001) 
suggest that > 20% of the 850 /xm emission in the 
outflow region could originate from line contami- 
nation. 

By contrast with the other two regions, the 
HH 24-26 (LBS 23) cores have relatively little 
extended submillimeter emission and are mostly 
compact (Mitchell et al. 2001; Johnstone et al. 
2001; Launhardt et al. 1996; Lis, Menten & Zylka 
1999) and cold (<10K, Chini et al. (1993)). Of the 
twelve condensations identified by Lis, Menten & 
Zylka (1999), our polarimetric image covers LMZ 
2, 3, and 4. All three of these cores have 3.6 
cm continuum sources (Bontemps, Andre & Ward- 
Thompson 1995; Gibb 1999). LMZ 3 (also known 
as HH24MMS, Chini et al. (1993)) is a class 
protostar (Bontemps, Andre & Ward- Thompson 



1995) , while LMZ 4 is a T Tauri star with a 
known CO outfiow (Snell & Edwards 1982; Gibb 
fc Heaton 1993). These are all indicators of low- 
mass star formation. 

NGC 2024 has an associated H II region and 
is the most prominent star formation region in 
Orion B, associated with a massive cluster, ion- 
izing B stars, and stars at all phases of evolution 
(Mezger et al. 1988; Lada et al. 1991b; Chandler 
& Carlstrom 1996). The submillimeter continuum 
emission was discovered by Mezger et al. (1988). 
The emission arises from a dense ridge of gas and 
dust behind the H II region, determined from the 
velocity of associated gas (Crutcher et al. 1986; 
Barnes et al. 1989), and consists of at least seven 
sources aligned along a ridge, similar to OMC- 
3 in Orion A (Johnstone & Bally 1999; Chini et 
al. 1997). Relatively more mass exists in the fil- 
amentary gas in the latter region. Two of these 
cores (FIR 4 and FIR 5) are the origins of unipo- 
lar molecular outflows, one of which is very highly 
collimated and very extended (Sanders & Willner 
1985; Richer et al. 1992; Chandler & Carlstrom 

1996) , while the FIR 6 core exhibits a compact 
outflow (Chandler & Carlstrom 1996) and con- 
tains a water maser (Genzel & Downes 1977), a 
signature of intermediate-mass protostars. The 
rest of the cores show no sign of star formation 
activity (Visser et al. 1998). Fukuda & Hanawa 
(2000) present a numerical simulation of triggered 
star formation along a filament by the expansion 
of an H II region and apply their analysis to NGC 
2024 in particular. In their picture, compression 
triggers two cores (i.e., FIR 4 and FIR 5) to col- 
lapse, and then at a later time, further collapse is 
triggered further up and down the filament. This 
sequence^ was observed by Chandler & Carlstrom 
(1996) in their study of outflows from NGC 2024, 
where the dynamical ages of outflows from FIR 5, 
FIR 4 and FIR 6 are 1.4 x 10* yr, 2.6 x 10^ yr, 
and 400-3300 yr. 

Magnetic fields play a crucial role in the process 
of star formation, through the magnetic support 
of molecular clouds, dissipation of angular momen- 
tum in accretion disks, and the generation of jets 
and outfiows (see Heiles et al. (1993) and refer- 
ences therein) . Polarized thermal emission at sub- 
millimeter wavelengths from aligned dust grains 
traces the direction of the magnetic field structure 
projected onto the plane of the sky (Hildebrand 
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1988). Absorption and scattering arc usually neg- 
ligible in the submillimeter. This is in contrast to 
optical and near infra-red polarimetry, where both 
of these processes can cause contamination. With 
the recent development of focal plane bolometer 
arrays equipped with polarimeters, sensitive imag- 
ing polarimetry in the submillimeter is now possi- 
ble. 

Polarimeters functioning at 100 /im (aboard the 
Kuiper Airborne Observatory), 350 /im (at the 
Caltech Submillimeter Observatory), and 850 /xm 
(at the James Clerk Maxwell Telescope) have de- 
tected polarized emission from dust toward many 
Galactic molecular clouds. These include the well- 
studied OMC-1 core in Orion (Schleuning 1998; 
Schleiming et al. 1997; Coppin ct al. 2000; Dotson 
et al. 2000); Sagittarius A (Aitken et al. 2000); 
Sagittarius B2 (Dowell et al. 1998) and recently 
OMC-3 in Orion (Matthews & Wilson (2000, here- 
after Paper I); Matthews et al. (2001, hereafter 
Paper II); Dowell 2001, in preparation). Nu- 
merous protostellar and starless cores have also 
been mapped (i.e., Holland et al. (1996); Ward- 
Thompson et al. (2000)). 

Polarimetry provides information only on the 
plane-of-sky magnetic field orientation, but no di- 
rect information about the magnetic field strength, 
since the degree of polarization is also dependent 
on other factors such as grain shape, degree of 
alignment, and composition. To obtain informa- 
tion about the magnetic field strength requires 
observation of Zeeman splitting of molecular or 
atomic spectral lines, which additionally provides 
information about the direction of the field along 
the line of sight. 

In this paper, we present the first submillime- 
ter polarimetry of the NGC 2071 core and the LBS 
23N region. Zeeman observations do not exist to- 
ward either of these two clouds. Far-infrared po- 
larimetry at 100 nm (Dotson et al. 2000) and 
Zeeman splitting observations of OH (Crutchcr & 
Kazes 1983; Kazes & Crutcher 1986; Crutcher et 
al. 1999) exist for the NGC 2024 ridge of cores. 
In § 2, we describe the 850 /im observations and 
the data reduction. The data are presented in § 3 
and the polarization patterns are interpreted in § 
4. We summarize our results in § 5. 



2. Observations and Data Reduction 

The observations presented here were obtained 
on 1998 September 7 and 8 using SCUBA (Sub- 
millimeter Common User Bolometer Array) (Hol- 
land et al. 1999) at the James Clerk Maxwell 
Telescope^ combined with the SCUBA polarime- 
ter (Greaves et al. 2000, 2001b). The nights were 
very stable, with t(225GHz) ranging from 0.05 
to 0.07 during the period of observations. Cal- 
ibration of the polarizer was performed on 1998 
September 5 using the Crab Nebula, for which a 
percentage polarization p = 19.3±4% and position 
angle ^ = 155±5° were measured. 

Observations were made of two overlapping 
fields toward NGC 2024, and one field each to- 
ward NGC 2071 and LBS 23N. Each polarization 
cycle consists of 16 integrations at 22.5° rotation 
intervals in 8 minutes integration time. Six obser- 
vations (48 minutes integration time) were made 
toward NGC 2071 and each position toward NGC 
2024, and 27 observations (3.6 hours integration 
time) were made toward LBS 23N. 

For each polarization cycle map, the raw 
SCUBA data were reduced using standard SCUBA 
software (Holland et al. 1999) to perform nod 
compensation, flatfielding and extinction correc- 
tions. More information on the polarimeter and its 
data cycle can be found in Greaves et al. (2000). 
For polarimetry data, extinction is estimated by 
extrapolation from the current CSO tau value at 
225 GHz. Bolometers with anomalously high noise 
were flagged, and the data were clipped at the 10 
(J level. Sky noise and instrumental polarization 
(IP) removal was then performed. The removal 
of sky noise is, even for basic SCUBA data, a 
task fraught with danger, but is even more so for 
polarimetry. The techniques of sky removal are 
discussed in detail in Paper II; we have taken care 
to select empty bolometers for sky subtraction. 

To first order, the removal of sky effects is 
achieved by chopping the JCMT secondary during 
SCUBA observations. Systematic errors can be in- 
troduced if significant polarized flux exists at the 
chop, or reference, positions. Large chop throws 

■^Thc JCMT is operated by the Royal Observatory Edin- 
burgh on behalf of the Particle Physics and Astronomy 
Research Council of the United Kingdom, the Netherlands 
Organization for Scientific Research, and the National Re- 
search Council of Canada. 
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of 150" were used for these observations (180" is 
the maximum for the JCMT) to try to chop as far 
off the bright emission as possible. The possible 
effects of polarized flux in the reference beams (or 
bolometers used for sky removal) are discussed in 
Paper 11. 

Since we have access to large-scale 850 /xm scan 
maps of Orion B north (Mitchell et al. 2001) and 
south obtained by the Canadian Consortium for 
Star Formation Studies at the JCMT, compar- 
isons can be made between intensities of candidate 
bolometers for sky noise subtraction in the on- 
source position and the fluxes in the corresponding 
chop positions. For two of our regions, NGC 2024 
and NGC 2071 , comparison reveals that the differ- 
ences between these two positions are significant, 
and that no bolometers on-source were "empty". 
The region around LBS 23N, on the other hand, 
reveals the same flux levels at the chop positions as 
at the edges of the on-source pointing. Thus, some 
bolometers at the edges of that map had average 
intensities of zero. To deal with sky bolometers 
which do not contain zero flux, the mean value 
per bolometer removed during sky noise subtrac- 
tion is added back into the maps of NGC 2024 
and NGC 2071, increasing the intensity in each 
bolometer. The main consequence of this addition 
is to reduce the polarization percentage, which is 
inversely proportional to /. 

After correcting for source rotation across the 
array, the Stokes' parameters were calculated by 
comparing measurements offset by 45° in wave- 
plate rotation (90° on the sky). The /, Q, and 
U maps for each pointing center were averaged, 
and standard dc^viations were derived by compar- 
ing the individual data sets. The two overlapping 
fields toward NGC 2024 were then combined into a 
mosaic; polarization data are uncalibratod, so the 
data were combined simply by averaging. Since no 
absolute flux calibrations are done on polarimetry 
data, the zero points are unknown. This presents 
problems for creating maps requiring more than 
one SCUBA field, since one does not know which 
field contains the best base level to use as refer- 
ence. We want to obtain the best estimates of 
uncalibrated flux without giving undue weight to 
either of our mapped fields. This was done by 
averaging the baseline levels in the two maps. 

The maps were then binned spatially by a factor 
of 3 (NGC 2071) or 4 (NGC 2024 and LBS 23N) 



in both RA and DEC to yield 9 or 12" sampling. 
Selection of "good" polarization vectors is done 
by filtering out all vectors for which p < 1% (since 
the IP values are accurate to just ±0.5%) and then 
filtering by signal-to-noise of polarization percent- 
age, uncertainty of polarization percentage, and 
the intensity level as appropriate to each field. 
The thresholding of p ensures that the effects of 
sidelobe polarization are low (see § 2.1). As fur- 
ther confirmation of the observed polarization pat- 
terns, we divided the data sets for each source into 
two subsets and performed complete reduction on 
each. The polarization patterns were qualitatively 
consistent with one another in all cases. 

We have truncated our selection of vectors at 
levels of 2% of the NGC 2071IR peak intensity, 
2.5% of the NGC 2024 source peak and 4% of 
the LBS 23N peak. Paper II illustrates that per- 
centage polarizations and position angles are reli- 
able down to these levels in cases where the ref- 
erence position's intensity is approximately 2% of 
the source peak intensity and the polarization per- 
centage in the reference position is no greater than 
that of the source field. Thus, in the lowest flux re- 
gions, the polarization percentage could be overes- 
timated by up to a factor of two while the polariza- 
tion position angle is incorrect by < ±10°. Since 
JCMT data obtained by chopping off the source 
record only the difference between these two posi- 
tions, our data provide no information about the 
polarized emission at the reference position. How- 
ever, based on more extensive intensity maps pro- 
duced by the scan-mapping technique, the chop 
positions are approximately 1%, 2% and < 1% of 
the on-source peaks for NGC 20711R, NGC 2024, 
and LBS 23N respectively. These are effectively 
lower limits to the fractional fluxes relative to the 
peaks, since scan mapping data may remove the 
fluxes of very extended features. If there is a uni- 
form background chopped out of the scan maps as 
well, then the reference position could in fact con- 
tain a more significant fraction of the source peak 
fiux. 

2.1. Sidelobe Polarization 

Even though the sidelobes of the JCMT beam 
contain less than 1% of the main beam power at 
850 /xm, significant polarization can be measured 
there because the two incoming planes of polar- 
ization experience different optical effects. This 
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can produce significant sidelobe polarization, that 
can become a source of systematic error especially 
when mapping extended fields in which the source 
located in the main beam is fainter than sources 
elsewhere in the SCUBA field. One can estimate 
the minimum believable source polarization per- 
centage from the expression: 

Pc..>2xp,(^)(§i) (1) 

where psi is the IP in the relevant part of the side- 
lobe, Psi /Pmb is the ratio of the sidelobe power to 
that in the main beam, and Fgi/F is the ratio of 
the mapped source flux in the sidelobe to that in 
the main beam (Greaves et al. 2001b). Basically, 
Pcrit is an estimate of the polarization percentage 
produced by the sidelobe source alone; hence only 
polarizations in excess of that value are believable. 

For the regions of NGC 2071 and LBS 23N, 
sidelobe polarization will be minimal, since in both 
those fields, the field centers were located at the 
submillimeter peaks. However, in NGC 2024, our 
northern mapping field was centered roughly on 
FIR 3, but FIR 5 was at the edge of the SCUBA 
map, approximately 70" from the center. The ra- 
tio of the FIR 5 flux to the central flux is ~ 3.4. 
Using a map of Saturn from 31 August 1998, the 
average power at 70" from the field center is ap- 
proximately 0.3% of the main beam. The aver- 
age instrumental polarization is 5.5%. Substitu- 
tion in equation (1) yields Pcrit > 0.11%. The 
planetary map was generated with a chop throw 
of 120", somewhat less than the 150" chop used 
in NGC 2024. However, using planetary data on 
Saturn from 11 October 1999 with a 150" chop, 
we obtain an estimate of pcrit > 0.2%. Hence, we 
are confident that our thresholding of polarization 
percentage at 1% removes all effects of sidelobe 
polarization in our maps. 

3. The Polarization Data 

The polarized emission in the submillimeter is 
thermal emission originating from rapidly spin- 
ning non-spherical grains, which are expected to 
be aligned, on average, perpendicular to the mag- 
netic field. The precise details of the alignment 
processes are not entirely understood. However, 
it is well understood that supra-thermally rotat- 
ing grains align perpendicular to the magnetic field 



on very short time-scales. Thus, grain alignment 
depends critically on the mechanisms by which 
grains are accelerated to supra-thermal rotational 
velocities. Radiative torques from short wave- 
length radiation (Draine & Weingartner 1996) and 
the well-known Purcell "rocket" mechanism (Pur- 
cell 1979) have both been proposed as candidate 
mechanisms. At present, however, these mecha- 
nisms appear insufficient to accelerate grains to 
the required rotational speeds in regions of high 
optical depth within molecular clouds (Lazarian 
et al. 1997). However, the data presented in this 
paper, as well as recent papers by other authors 
(Paper II; Dotson et al. (2000); Schleuning (1998)) 
provide evidence that grain alignment does oc- 
cur in dense regions of molecular clouds and that 
the underlying field structure is well-ordered, as 
demonstrated by the highly structured polarized 
emission seen in our maps. Johnstone et al. (2001) 
compare 850 /zm data of NGC 2071 and LBS 23 
to CS J = 2 — 1 maps by Lada et al. (1991a) and 
conclude that these clumps of dust emission occur 
within regions above a column density threshold 
of N{H) ~ 10^^ cm~^, which corresponds to an 
Ay > 5. Thus, the 850 /im dust emission is asso- 
ciated with extinctions exceeding those at which 
radiative alignment is expected to be significant. 
This result poses an important challenge for grain 
alignment models. 

The shapes of grains, their size distribution, 
their composition, and the details of their align- 
ment mechanisms all strongly affect the ability 
of grains to produce polarized thermal radiation, 
and hence influence the resulting polarization pat- 
terns. None of these factors are well-constrained, 
either theoretically or observationally. However, 
Fiege & Pudritz (2000c) developed a pragmatic 
approach that combines all of the grain and align- 
ment properties for a distribution of grain species 
into a single parameter (a), which can be ob- 
tained observationally. This parameterization per- 
mits modeling of polarization patterns without de- 
tailed knowledge of grain and alignment physics. 
It was assumed in Fiege &: Pudritz (2000c) that all 
grain species are uniformly mixed with the gas and 
aligned to the same extent everywhere within the 
cloud, which implies that (a) is constant within 
any given cloud. However, their model is easily 
extended to allow systematic variations in grain 
and alignment properties throughout the cloud, 
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by simply allowing (a) to be a function of ei- 
ther the density, or the distance from the centre 
of the cloud. We find that allowing (a) to vary 
within a cloud docs not typically have a large im- 
pact on the geometry of the polarization patterns, 
although the polarization percentage and the de- 
gree of depolarization toward bright sources can 
vary more substantially. In this sense, our models 
are quite robust against changes to the underlying 
grain composition and alignment properties. 

If we assume, for now, that (a) is constant, 
then the maximum possible polarization percent- 
age (for the most favourable geometry) with a con- 
stant field parallel to the plane of the sky, is given 
by the formula 

^™""~l-(a>/6 ^' 

(Fiege & Pudritz 2000c) . The maximum polariza- 
tion percentages for the three regions that we have 
mapped are 15% for NGC 2071, 13% for NGC 
2024, and 22% for LBS 23N. The value for LBS 
23N may be anomalously high; the next highest 
percentage measured is 12.6%. Therefore, 15% is 
a reasonable maximum value for NGC 2071, with 
13% an upper limit for NGC 2024 and LBS 23N. 
Using 15% as an upper limit on the polarization 
percentages in these regions, and assuming this is 
close to the theoretical maximum for a line of sight 
passing through a magnetic field with ideal geom- 
etry, then we find that (a) w 0.14 for all three 
regions. This estimate is really a lower bound for 
(a) because the observed polarization percentage 
is reduced by any component of the field lying 
along the line of sight. Nevertheless, we are en- 
couraged that our estimates of (a) do not vary 
enormously between the three regions, and con- 
clude that a = 0.14 is a reasonable global esti- 
mate for (a) in Orion B. We use this value of (a) 
to generate models of NGC 2024 (see §4.3). 

3.1. Polarization Patterns 

3.1.1. NGC 2071 

In Figure 1, we present the 850 /um polariza- 
tion vectors of NGC 2071 IR, superimposed on 
a false color image of the 850 /xm dust emission 
(from Mitchell et al. (2001)). The inset shows a 
"blow-up" of the central region in order to show 
the weak polarization features better. The data 



are also presented in Table 1 of Appendix A. 

Toward the region of strongest emission (the 
IR cluster), the percentage polarization is quite 
weak, typically less than 2%. However, polariza- 
tion as strong as 15% is measured toward fainter 
parts of this region. The polarization pattern is 
highly symmetric, including the curvature of the 
vector orientation around the peak to the north- 
east and southwest. Over the entire core region, 
the vectors exhibit a mean orientation of 20° east 
of north. Based on low-flow velocities of HCO+ 
maps of the dominant outflow in the region, Gi- 
rart et al. (1999) estimate its overall position an- 
gle to be » 40° (east of north). Figure 1 shows 
this outflow orientation schematically. Closer to 
the emission peak, the dominant position angle is 
~ 30°, which is within 10° of the outflow's orien- 
tation. The outflow is thus well-aligned with the 
polarization vectors. 

Qualitatively, the polarization pattern is similar 
to that of OMC-1, which was interpreted as evi- 
dence for an hourglass magnetic field geometry in 
that region (Schleuning 1998). In this picture, the 
magnetic field is pinched toward the central source 
due to flux freezing with the infalling gas. How- 
ever, we do not observe the same sort of flattening 
of the NGC 2071 core that is evident in the 350 
lj,m flux map of OMC-1 (Schleuning 1998). The 
NGC 2071 core maintains a circular shape even far 
from the peak, with no suggestion of oblateness or 
prolateness (see Figure 1). 

3.1.2. LBS 23N 

Figure 2 shows the dust polarization of the LBS 
23N source and its surrounding region (the data 
are also presented in Table 2 of Appendix A). A 
string of cores is surrounded by a faint filamentary 
background of dust emission. The dominant orien- 
tation of the polarization pattern is north-south, 
except at the southern boundary, where the vec- 
tors lie almost east- west. The strongest polariza- 
tion (22%) is seen in the east and west limbs of 
the "filament". 

It is intriguing that the vectors abruptly change 
direction by 90° near the southern- most boundary 
of our polarization map to become aligned east- 
west. The 850 /xm continuum map of Mitchell et 
al. (2001) shows that the underlying dust struc- 
ture widens south of our polarization map to a fil- 
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amcntary segment oriented roughly east-west (see 
Figure 2). A similar shift of polarization orienta- 
tion has been observed in the southern part of the 
OMC-3 filament in Orion A, where the filament's 
projected width on the plane of the sky appears to 
increase (Papers I and II). We discuss three possi- 
ble interpretations of our data in § 4.2. Polarimo- 
try in the southern region of LBS 23 will help us 
in our modeling effort to understand the struc- 
ture of the magnetic field in this region. However, 
the faintness of the LBS 23 filament will make ob- 
taining the necessary data more difficult than in 
OMC-3. 

3.1.3. NGC 2024 

The polarization pattern of the molecular ridge 

of NGC 2024 is shown in Figure 3, superimposed 
on a false color image of the 850 /xm dust emis- 
sion.^ The data are presented in Table 3 of Ap- 
pendix A. Seven continuum sources (FIR 1-7) 
have been identified along the ridge (Mezger et 
al. 1988, 1992). In the north (FIR 1-4), the 
vectors along the peaks arc weak and parallel to 
the ridge. The southern sources (FIR 5-7) have 
very weak polarization toward the brightest cores 
and are significantly depolarized. To the east and 
west, the polarization becomes quite intense (up 
to 10%), and a change in vector orientation of ap- 
proximately 90° is observed from west to east. 

Linear polarization at 100 /um has been ob- 
served previously toward this region with the 
Kuiper Airborne Observatory (Hildebrand et 
al. 1995; Dotson et al. 2000). These data are 
centered on the FIR 5 core and cover the same 
area as our map. Figure 4 shows our polarization 
data and those of Dotson et al. (2000) at 100 /zm 
plotted over contours of the 850 nm total emission. 
The 100 /im data are more sparsely sampled than 
the SCUBA data (owing to the 35" beam), but the 
two patterns are remarkably consistent. The vec- 
tors trace the same V-shaped pattern seen at 850 
/um, although the finer sampling by SCUBA makes 
the pattern much easier to see. The chop throw of 
the KAO observations was 10.5', over four times 
greater than that used for our data set. The fact 
that the polarization patterns are very consistent 
suggests that the position angle at least has not 

^JCMT data taken by Canadian Consortium for Star- 
Formation Studies, as yet unpublished. 



been strongly affected by any systematic effects 
due to chopping. The maximum polarization per- 
centage measured at 100 /xm is 5.86l5'oo%j while 
the minimum vector has a polarization percentage 
of 0.36t[J;g^% and a position angle of 164.9 ± 5.2° 
(east of north), where the uncertainties are abso- 
lute upper limits. At roughly the same position, 
the 850 /zm polarization is 1.63±0.15% at an angle 
of -20.7 ± 2.7° (equivalent to 159.3 ± 2.7°). This 
is a relatively low polarization percentage at 850 
/xm, and the position angles of the vectors agree 
within measurement uncertainties. At 850 /xm, 
the maximum and minimum polarization percent- 
ages are 13.41 ± 1.13% and 1.02 ± 0.09%. Where 
the northern and southern parts of the filaments 
meet (i.e., between FIR 4 and 5), there is signifi- 
cant depolarization regardless of the low fiux level. 
Depolarization in the filamentary material away 
from the cores is also observed between cores in 
OMC-3 in Orion A (see Paper II). 

There is no obvious correspondence between the 
polarization position angles and the directions of 
outfiows in the region. The strongest outflow is 
unipolar and extends from FIR 5 to the southeast 
(Richer, Hills & Padman 1992). Greaves, Holland 
& Ward- Thompson (2001) have measured polar- 
ized spectral line emission from the CO J = 2 — 1 
line utilizing the Goldreich-Kylafis effect (Goldre- 
ich & Kylafis 1981) in the outflow of FIR 5. The 
polarization vectors measured in the outflow are 
plotted on Figure 3, along with the orientation of 
the outflow itself. For large CO optical depths, 
line polarization is parallel to the fleld, unless the 
field is at a large angle to the velocity gradient 
(Kylafis 1983). Therefore, Greaves et al. (2001a) 
conclude that the magnetic field associated with 
the outflowing gas is aligned with the outflow it- 
self. If the large scale fleld in the region main- 
tains the same orientations, we would expect the 
continuum polarization vectors to lie perpendicu- 
lar to the outflow. But the continuum vectors do 
not exhibit this behavior. Nor does their orien- 
tation appear to be influenced by their proximity 
to the outflow. The dust polarization data do not 
exhibit an obvious correlation with either the out- 
flow or the line polarization data which traces its 
magnetic fleld. Thus, we conclude that a complex 
field geometry threads this region. 
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3.2. Depolarization toward High Intensi- 
ties 

As observed in many other regions, a declining 
polarization percentage is detected toward regions 
of increasing intcuisity in each of our fields. Figure 
5 shows the percentage polarization as a function 
of unpolarized I for each region. The depolariza- 
tion effect is clearly a global feature in our maps 
and is not limited to a single region or core. The 
slopes of Figure 5 range from —0.5 to —0.9, but 
the trend of declining polarization percentage is 
clear. The same trend was measured in the OMC- 
3 region of Orion A (see Paper II). It is possible 
for a "polarization hole" to result from a system- 
atic error produced by chopping to remove back- 
ground sky during observing. However, Paper II 
shows that depolarizations as steep as those mea- 
sured in Orion B (i.e., slopes of log / versus log 
p < —0.55) cannot be produced by the effects of 
chopping alone, particularly at high values of /, 
even for the case where the flux at the reference 
position is 25% that of the source peak and the 
reference polarization percentage is twice that of 
the source. Therefore, there must be a physical 
cause for this effect. 

The issue of depolarization in dense regions 
is particularly interesting as a test of alignment 
mechanisms. Grain alignment theories at present 
cannot explain why grains arc so well aligned in 
regions of high optical depth within molecular 
clouds. For example, the Purcell mechanism (Pur- 
cell 1979) and its variations are not likely to work 
in these regions, where atomic hydrogen is almost 
completely absent (Lazarian et al. 1997). The 
radiative alignment mechanism (Draine & Wein- 
gartner 1996), which aligns grains efficiently in the 
ISM, is unlikely to work in dense molecular re- 
gions, where short wavelength starlight (needed to 
spin up and align grains) is extinguished. Appeal- 
ing to embedded sources to provide the necessary 
radiation is not likely to help, since this would pre- 
dict the best alignment to be associated with the 
dense regions around evolved protostars, while the 
opposite is observed. Thus, our observations pose 
an interesting and extremely important challenge 
for grain alignment theory. 

We consider 3 possible explanations for the "po- 
larization hole" effect. First of all, it is possible 
that the depolarization in our maps are the result 



of grains in the central regions of highest density 
that are either poorly aligned, or intrinsically poor 
polarizers (Weintraub, Goodman fc Akeson 2000; 
Goodman et al. 1995). In such a case, all dust 
grains through the cloud contribute to the con- 
tinuum intensity, but only the fraction that are 
aligned contribute to the polarized flux given by 
the Stokes' vectors Q and U. Therefore, toward 
regions of high intensity, the ratio of polarized to 
total intensities is systematically reduced. Alter- 
natively, Fiege & Pudritz (2000c) proposed that 
the depolarization could result from the large scale 
field structure itself, since helically twisted fields 
result in depolarization as a result of oppositely 
signed contributions to the Stokes' vectors from 
the poloidal and toroidal field components. This 
field geometry has been used to model the OMC-3 
part of Orion A's Integral-shaped Filament (Paper 
II), and we show a similar model for NGC 2024 in 
§ 4.3.1 below. A third possibility is that a dis- 
ordered, or tangled, field exists on scales smaller 
than our beam, and the field averages to zero in 
the JCMT's 14" beam. It is possible that more 
than one of these effects could be contributing to 
the depolarization observed in our maps. 

Interestingly, NGC 2071 does not exhibit a po- 
larization hole coincident with the peak of the dust 
emission; instead, two polarization holes are de- 
tected, offset to the northeast and southwest of the 
dust emission peak. These polarization holes are 
not due to bad bolometers or low signal-to-noise. 
The degree of polarization in those locations is less 
than 1% and is thus considered negligible at these 
positions. In NGC 2024, depolarization is so sig- 
nificant in the case of FIR 5-7 that polarization is 
u(-gligible across those cores. In LBS 23N, depolar- 
ization can be seen in Figure 2 toward the bright- 
est core, OriBN 59 (see Mitcheh et al. (2001)). 
It is difficult to sample the depolarization in the 
fainter sources of Figure 2 because they are mire- 
solved, but OriBN 58 does show significant depo- 
larization. Our mapping does not extend to the 
OriBN 60 core, and further data toward LBS 23 
(south) will be needed to determine whether or not 
the observed depolarization extends to this region. 

4. Interpreting Polarization Patterns 

Polarization vectors in the submillimeter are 
perpendicular to the direction of the magnetic 
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field component in the plane of the sky, appro- 
priately averaged along the line of sight through 
the cloud. However, it is generally insufficient to 
infer the field structure by simply rotating the po- 
larization vectors by 90°, since this procedure im- 
plicitly assumes that the direction of the plane of 
sky component of the magnetic field is constant. 
Magnetic fields in molecular clouds are almost cer- 
tainly curved in three dimensions, which violates 
this condition for all but the simplest geometries. 
A simple example serves to illustrate that this 
technique fails. Fiege & Pudritz (2000c) recently 
calculated the polarized emission from their mod- 
els of filamentary clouds threaded by helical mag- 
netic fields (Fiege & Pudritz 2000a). This was 
done by numerically integrating the contributions 
to the Stokes' vectors along lines of sight through 
the filament. They demonstrated that the po- 
larization is always parallel or perpendicular to 
the symmetry axis of the filaments in the plane 
of the sky, (sometimes with sudden 90° flips in 
orientation at some radial distance). Now, if we 
try to interpret the resulting polarization maps by 
simply rotating the polarization vectors derived 
from their models, we discover that this approach 
does not resemble the input helical field geome- 
try. A more sophisticated approach is generally 
required to compare candidate models, threaded 
by non-trivial 3-dimensional magnetic fields, with 
the data. This must involve direct modeling of 
polarization patterns. 

It is important to note that interpretations of 
polarization data are somewhat degenerate. By 
varying the magnetic field geometry and the po- 
larization efficiency of the grains, it is often possi- 
ble to produce a qualitatively similar polarization 
map from more than one field geometry. An ex- 
ample of this difficulty is given in § 4.3, where 
we present two alternate models of our NGC 2024 
polarization map. It is generally impossible to 
uniquely determine the 3-dimensional structure of 
the magnetic field from polarization data alone. 
Nevertheless, polarization maps can be used to 
provide strong constraints on models of magne- 
tized clouds, filaments, and cores as follows. 1) 
On the basis of polarization data alone, it is often 
possible to rule out models that are inconsistent 
with the data when viewed from all possible ori- 
entations. For example, in Paper II, we were able 
to rule out all models of the Integral-shaped Fil- 



ament in Orion A that are threaded by a purely 
poloidal field. 2) All existing magnetohydrostatic 
models of clouds, filaments, and cores are axisym- 
mctric. The assumption of axisymmetry reduces 
the set of models to be compared with the data 
which limits the allowed configurations. We find 
that much of the polarization degeneracy is lifted 
when only axisymmetric models, or perturbations 
of axisymmetric models (see § 4.3.1), are consid- 
ered. We also find that certain field structures 
produce identifiable, although sometimes surpris- 
ing, "signatures" in the polarization data, which 
can be used guide toward developing detailed 
models. 3) There are few regions where both po- 
larization data and Zeeman data exist. Zeeman 
data can be combined with the polarimetry to 
construct a more nearly unique model than could 
be determined by polarimetry alone, since Zee- 
man observations measure the line of sight compo- 
nent of the magnetic field. However, one must be 
extremely careful because Zeeman measurements 
and submillimeter polarimetry usually trace differ- 
ent components of the gas. An interpretation that 
combines Zeeman data and polarimetry implicitly 
assumes that the field is related in these compo- 
nents. We combine polarization and Zeeman data 
in § 4.3.2, where we model the NGC 2024 as a shell 
swept up by an H II region around a dense ridge. 

An alternative approach is to compare the data 
with polarization maps predicted from simula- 
tions. The cores arising in simulations are not 
generally in equilibium, and the field simulations 
are not generally axisymmetric. This technique 
has the potential to provide useful insight into the 
general statistical behavior of polarization maps 
for ensembles of objects. Padoan et al. (2001) 
have applied the polarization modeling technique 
developed by Fiege & Pudritz (2000d) to simu- 
lations of MHD turbulence. They find that the 
relationship between polarization percentage and 
submillimeter intensity is a robust statistic, which 
can be used to constrain models of grain align- 
ment. A limitation of turbulence simulations is 
that they cannot be easily used to model specific 
objects, which is the focus of this paper. 

4.1. NGC 2071 

NGC 2071 is a massive core that has no asso- 
ciated filamentary structure. Nevertheless, there 
are some strong similarities between this core and 
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OMC-1, even though the latter is associated with 
the massive Integral-shaped Filament of Orion A. 
There are no Zeeman data toward NGC 2071, so 
the polarization data alone constrain the field ge- 
ometry. The polarization map is similar to the 100 
/Ltm polarization map of OMC-1 made by Schleun- 
ing (1998). In that region, the polarization data 
are interpreted as indicating an "hourglass" mag- 
netic field, which developed as the core contracted 
perpendicular to field lines. The contraction is 
centered on the infra-red source known as KL. 
However, there is an important difference between 
Figure 1 and Schlcuning's 100 and 350 /im maps. 
The OMC-1 core exhibits a flattened structure, 
whose long axis (in projection) is parallel to the 
mean direction of the polarization vectors. This is 
precisely what one would expect for an oblate core 
threaded by an hourglass field (Mouschovias 1976; 
Tomisaka, Ikeuchi, & Nakamura 1988). However, 
the NGC 2071 core is not flattened at all. There 
is no evidence of elongation either parallel or per- 
pendicular to the outflow in our data (see Figure 
1) or the 50 and 100 /xm observations of Butner et 
al. (1990). 

The direction of the most powerful outflow of 
NGC 2071, thought to originate from the source 
IRS 3, is shown in Figure 1. It is aligned at a 
position angle of 40° (east of north). The mean 
polarization position angle in our data set is ~ 
22°, with a dispersion about the mean of 32°. In 
the subset of vectors toward the peak, the mean 
position angle is 34° with a standard deviation 
about the mean of 33°. The vectors are generally 
consistent with with the direction of the outflow. 

Houdc ct al. (2001) have recently used a new 
technique, which compares the line profiles of 
the coexisting neutral (HCN) and ionic (HCO+) 
species, to show that the outflow's magnetic field 
and its direction are well-aligned in NGC 2071. 
However, our polarization data suggest a magnetic 
field that is either orthogonal to the outflow direc- 
tion, or strongly toroidal about an axis parallel to 
the outflow. In either case, the magnetic field in 
the core does not appear to coincide with the field 
in the outflow. The first possibility would be diffi- 
cult to reconcile with theory, unless the magnetic 
field at the small scales where the outflow origi- 
nates is unrelated to the large scale fleld thread- 
ing the core. It would be especially problematic 
for models of protostellar formation that favor col- 



lapse along field lines, resulting in outflows aligned 
with the field of the core. This possibility would 
substantially diminish the role of magnetic fields 
in protostellar collapse. The latter possibility is 
more intriguing because it could be reconciled with 
a class of self-similar outflow models proposed by 
Fiege & Henriksen (1996) and elaborated on by 
Lery et al. (1999). These models permit a mag- 
netic fleld which is primarily toroidal in the dense, 
partially pressure supported, gas near the equa- 
torial plane, even far from the central star. The 
topology of the field is quadrupolar in the poloidal 
plane, so that the field is almost radial near the 
symmetry axis defined by the outflow direction 
and parallel to the outflow. 

Motte et al. (2001) have compared the SCUBA 
850 fim emission from NGC 2071 to pre-existing 
maps of CO J = 3 - 2 and HCO+ J = 4 - 3 
(Chernin & Masson 1992) and find that up to half 
the emission in the 850 /xm bandpass could be orig- 
inating in the outflow's line emission. This could 
affect the interpretations given above. Motte et 
al. (2001) estimate that > 20% of the contin- 
uum emission at 850 /xm could arise from the CO 
J = 3 — 2 line. According to theory, in the case of 
high density gas, the linear polarization produced 
in spectral lines should be parallel to the mag- 
netic field. Even if the magnetic fields giving rise 
to the dust and line polarizations were perfectly 
aligned, they would contribute orthogonal vectors 
to the net polarization. Thus, the usual interpre- 
tation of submillimeter continuum polarization as 
being orthogonal to the plane-of-sky component 
of the magnetic fleld could be incorrect in regions 
where the CO emission is strong enough to dom- 
inate the continuum flux from dust. If our polar- 
ization data were dominated by line emission, then 
the inferred field direction would align to within 
10° of the CO outflow. However, where the CO 
emission is optically thick, the maximum polar- 
ization contributed by the outflow is on the order 
of 1%, as is the case for CO J = 3 — 2 emission 
from the outflow of FIR 5 in NGC 2024 (Greaves et 
al. 2001a). Therefore, when dust polarizations ex- 
ceed 1%, they are unlikely to be dominated by CO 
line contamination (J. Greaves, private communi- 
cation). Therefore, we conclude that the strong 
polarizations measured toward this core are dom- 
inated by dust emission and hence do not infer a 
magnetic field aligned with the outflow direction. 
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Finally, we note that the orientation of the out- 
flow relative to the core's magnetic field may turn 
out to be inconsequential. NGC 2071 is a mas- 
sive core forming many stars. This is very differ- 
ent from the case of low mass cores forming iso- 
lated stars, to which most of the existing theory 
applies. Within several star-forming regions, mul- 
tiple outflows of various orientations are observed, 
i.e., OMC-2/3 in Orion A (Aso et al. 2000; Yu et 
al. 2000) . If the collapse of massive cores is really 
dominated by the mean magnetic field, then one 
might expect a rough alignment between all out- 
flows originating from the same core. Since this 
is not observed in NGC 2071 (Eisloffel 2000), it 
is unclear how relevant the large-scale field is in 
determining the orientation of outflows. 

4.1.1. An Estimate of Magnetic Field Strength 

Chandrasekhar & Fermi (1953a) demonstrated 

that when a mean field direction can be defined, 
the magnetic field strength can be derived using 
the dispersion in polarization angles, ag, the dis- 
persion in velocity along the line of sight, cr^,,^^, 
and the mean density, p, in the region local to the 
polarization and velocity data. In cases of com- 
plex, ordered field geometries, it is no easy task to 
assign a mean field direction based on polarization 
data (see §§4.2 and 4.3 below). However, in NGC 
2071, we can apply the method of Chandrasekhar 
& Fermi (1953a) assuming that the field thread- 
ing the core is relatively straight. The mean field 
strength, {B), is given by 
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(3) 



The NGC 2071 core was mapped in NH3 (1,1) 
emission at high velocity resolution by Takano 
et al. (1986). The ammonia emission extends 
over 2' from the dust peak location and shows 
more elongation than the dust map. (Takano ct 
al. 1986) identified three components in the line 
profiles from across the core. The narrowest fea- 
ture is described as a "spike" in the profile and is 
seen over the entire mapped area. This emission 
is attributed to the core itself (as opposed to any 
outflowing gas). The FWHM of this component is 
estimated to be 1 km s~^, which corresponds to a 
dispersion of approximately 0.5 km s~^. Ammo- 
nia traces relatively high density gas (n(i?2) > 10* 



cm~^), so wc use 10 cm~^ as an estimate of the 
mean density across the core. As discussed above, 
the dispersion in polarization angles at 850 /itm 
is 33°, which corresponds to 0.58 radians. Sub- 
stitution of these values into equation (3) yields 
(B) » 56/LtG which is consistent with what would 
be expected within a dense core. 

Heitsch et al. (2001) propose an amendment to 
the Chandrasekhar and Fermi model which yields 

an estimate of the rms field strength. In the case 
of equipartition between the magnetic and kinetic 
energies, the total magnetic energy is given by 



[l + (3a,^)] 



(4) 



which gives similar results to equation (3) when 
the field is highly ordered. For NGC 2071, wc 
find that the ratio of the total magnetic energy, 
{B^), to the mean field energy, (S)^, is 2. Using 
the values for NGC 2071 above, the rms field is 
estimated to be 78 p-G, a factor of 1.4 times greater 
than the mean field. 

Crutcher (1999) compared the densities within 
star-forming clouds with the line-of-sight field 
strengths measured using the Zeeman effect on 
molecular lines. Based on his compilation of all 
existing data, he concluded there was no substan- 
tial evidence against equipartition in star-forming 
clouds. The line-of-sight field strength. Bios, can 
be extrapolated from the mean field strength 
based on a statistical average over all potential 
inclinations of the field to the line of sight. We 
estimate Bios to be roughly the mean field divided 
by 2, or 28 pG. Comparing this field strength to 
the density within the NGC 2071 core of 10"* cm'^ 
using Table 1 of Crutcher (1999), it is clear that 
this value is comparable to measurements of the 
Barnard 1 dark cloud core, but substantially less 
than the more massive M17 cloud. 

Based on this comparison, we conclude that, 
while the Chandrasekhar-Fermi method does give 
a value of mean field strength consistent with those 
measured directly using Zeeman splitting, it also 
predicts that approximately half the magnetic en- 
ergy does not reside in the mean field component. 
These results are not surprising given the fact that 
an ordered pattern is observed in the polarization 
data which is not straight. Thus, we should not be 
too quick to attribute the rest of the magnetic en- 
ergy to turbulence. However, direct measurement 
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of Zeeman splitting of OH lines at several posi- 
tions in NGC 2071 would provide a further test 
on the validity of a simple, mean field geometry 
within this cloud. 

4.2. LBS 23N 

The polarization structure of the LBS 23N re- 
gion is somewhat less orderly than that of either 
NGC 2071 or NGC 2024, and there are no Zeeman 
data toward this region. However, there remains 
sufficient structure in the polarization map to war- 
rant some comparison to filamentary models. The 
LBS 23N region is a linear string of roughly equally 
spaced condensations, which might be the result 
of fragmentation due to gravitational instabilities 
(see Chandrasekhar & Fermi (1953b), Nakamura, 
Hanawa, & Nakano (1993) and Fiegc & Pudritz 
(2000b) for a discussion of gravitational instabil- 
ities in magnetized filaments). The overall ori- 
entation of the polarization vectors is along the 
filament, with some degree of scatter in position 
angle. The mean position angle in the 70 vectors 
is 8.7° with a standard deviation of 5.0°. 

We offer three possible explanations for such a 
pattern. The first possibility is that the original 
filament is threaded by a magnetic field transverse 
to the symmetry axis. This magnetic field geom- 
etry would reproduce the orientation of the ma- 
jority of the polarization vectors. However, it is 
not easy to understand how such a field geome- 
try could be associated with a filament or linear 
string of cores collapsed in three dimensions, since 
a transverse field would more likely be associated 
with a sheet-like distribution of material. This 
geometry also cannot readily explain the 90° flip 
seen in the polarization vectors toward the south- 
ern end of the filament. 

Alternatively, Fiege & Pudritz (2000c) pro- 
posed that such a pattern could be explained 
as the result of a helical field dominated by the 
toroidal field component (Type 1, according to 
their classification) as observed in the OMC-3 fil- 
ament of Orion A. Helical field geometry can cre- 
ate polarization patterns where vectors flip by 90° 
at the interface between poloidally and toroidally 
dominated regions, which could occur if the field 
threading the southern part of the filament is less 
tightly wound than in the north. Such a flip is 
suggested by our data at the southern end of the 
polarimetry map. We note that these data pro- 



vide only sparse coverage of this region and are 
the lowest signal-to-noise data of all three regions. 
Modeling of this region will require observations 
of the southern part of the filament. 

Finally, Paper H presents a model for two 
crossed filaments, unconnected but overlapping in 
projection on the plane of the sky. In such a case, 
the vectors in the overlap region can be dominated 
by one of the two filaments. The east-west fila- 
mentary structure is identified in Figure 2. As for 
OMC-3, this scenario can be tested by more ex- 
tensive mapping into LBS 23 south. 

4.3. NGC 2024 

The overall appearance of NGC 2024's submil- 
limeter emission is that of a short filament, which 
is broken into three main condensations. The 
polarization vectors show a complicated but ex- 
tremely well-ordered pattern throughout the re- 
gion. An interesting feature of the map is a rather 
abrupt change in orientation across the central 
ridge of the filament. This feature cannot be rec- 
onciled easily with a purely poloidal field along the 
filament axis, since that geometry would produce 
polarization vectors perpendicular to the axis, ir- 
respective of the inclination of the filament. Fiege 
& Pudritz (2000c) found that polarization vectors 
can flip by 90° across the central axis of a filament 
threaded by a helical field, when the filament is in- 
clined at an angle relative to the plane of the sky. 
However, the polarization vectors in their model 
are always aligned either parallel to or perpendic- 
ular to the filament. This is not the case in our 
map, where the vectors form a "V" like pattern, 
with the apex located on the central ridge of the 
filament. However, we show below that a simple 
extension to their model can account for the qual- 
itative features of this map. 

The interpretation of the polarization pattern is 
constrained by measurements of the line-of-sight 
magnetic field component. Bios (sometimes de- 
noted i3||). Crutcher et al. (1999), using Zee- 
man splitting of OH absorption lines, detected a 
peak in Bios of almost 100 iiG at a position west 
of FIR 4. The field strength declines smoothly 
from this southwestern position in their map to 
the northeast, where it drops to zero (cf. their 
figure 3). OH samples gas densities of approxi- 
mately 10^ cm~^, which are substantially lower 
than the volume densities estimated along the 



12 



NGC 2024 molecular ridge. The extended fila- 
mentary structure in which the cores are embed- 
ded have < riH >~ 5 x 10® cm~^, while the cores 
themselves have < nn >~ 2 x 10*^ cm^^ (Mezger 
et al. 1992). Thus, it is not certain that the 
field being traced by the Zeeman data threads the 
same dense gas whose magnetic field is traced by 
our polarization measurements. However, the fact 
that the positions of low polarization percentage 
in both 100 /v,m and 850 /im correspond closely 
with the maximum in line-of-sight field could sup- 
port a correlation in the fields traced by the two 
techniques, since low polarizations could indicate 
a region where Bios 

is dominant. Conversely, the 
low polarization could be a produced by variations 
in alignment or grain properties and thus be un- 
related to Biog. Keeping this uncertainty in mind, 
we offer two theoretical models for this pattern. 

Schulz et al. (1991) measured NH3 and CS gas 
toward NGC 2024 and found distributions of dense 
gas that closely mirrored the pattern of dust con- 
densations. This supported later conclusions that 
the dense ridge is located on the far-side of the H II 
region (Chandler & Carlstrom 1996; Crutcher et 
al. 1999) ; as originally schematically illustrated by 
Barnes et al. (1989). We illustrate the most likely 
configuration in Figure 6. In order for OH to be 
observed in absorption, there must be a source of 
background continuum emission. Thus, either the 
OH Zeeman data trace the field in the foreground 
dust lanes (with the observable H II region as the 
continuum source) or the OH absorption is occur- 
ing in the denser molecular cloud, thus requiring 
a second source of continuum emission behind the 
ridge. For instance, the H II region could have 
penetrated to different depths of the cloud due to 
density variations. The latter scenario is favored 
by Crutcher et al. (1999) with the ridge located 
on the far side of the H II region. 

In the first case, the Zeeman splitting OH ab- 
sorption data may be sampling an envelope of gas 
either in the foreground of the H II region itself, 
or simply at the boundary between the H II region 
and the very dense molecular cloud behind it, re- 
gardless of the location of the dense ridge. In this 
case, the Zeeman data are unrelated to the polar- 
ization data and would not provide a constraint. 
In the second case, we assume that although the 
OH data sample a less dense region of gas, the 
measurements are sampling the same field geome- 



try in both regions. 

4.3.1. A Helical Field Model 

Perfectly straight, infinite filaments are an ide- 
alization, which may need to be relaxed somewhat 
when comparing models with real data. Figure 7 
shows a finite segment of a filament model from 
Fiege & Pudritz (2000a), which we have "bent" 
into a circular arc using a simple mathematical 
transformation that maps planes parallel to the 
axis of the filament into concentric cylinders of 
some specified mean radius of curvature, orthog- 
onal to the filament axis. The density and mag- 
netic field are transformed self-consistently, using 
the Lagrangian formulation of the MHD induction 
equation and the equation of mass conservation 
(see Parker for example): 

P Po 

where J is the Jacobian of the coordinate transfor- 
mation and |J| is its determinant. We will thor- 
oughly discuss the details of our modeling tech- 
nique in a future paper. 

The polarization structure of the theoretical 
map reproduces most of the qualitative features 
of Figure 3. On the other hand, models that are 
threaded by purely poloidal fields do not agree 
with the polarization data. Note that the model 
shown is a preliminary result and should not be 
regarded as our final, best fit to the data. Never- 
theless, we have provided the exact parameters of 
this model in the caption of Figure 7. 

It is interesting that a bend in the filament 
is required to reproduce the polarization pattern 
shown in Figure 3. If the filament is bent so that 
the source FIR 5 is closest to the expanding ion- 
ization front of the H II region, then this could ex- 
plain why it is the most evolved source along the 
ridge. The outflow from this source is the most 
powerful in the region (Richer et al. 1992), while 
evidence for less evolved outflows can be found in 
the sources FIR 4 and FIR 6 (Chandler & Carl- 
strom 1996). 
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4-3.2. Expansion of the H II Region around A 

Dense Filament 

The expected Zeeman pattern from a helical 
field is a reversal in orientation of Bios across the 
symmetry axis of a filament. If the field thread- 
ing the filament is related to the field in the lower 
density gas, then the helical field interpretation 
would not be consistent with the observations of 
Crutcher et al. (1999) toward the northern por- 
tion of the NGC 2024 dense ridge. In order to 
reconcile the Zeeman and polarization data, we 
suggest that the magnetic field has been distorted 
from its original structure by the expansion of the 
H II region, as shown in the schematic diagram of 
Figure 8. The field may be pushed out on the sur- 
face of the ionization front, but will bend around 
the dense, filamentary system of cores comprising 
the NGC 2024 region. 

The asymmetric Zeeman data requires one 
of the following geometries: either the ionizing 
source is in front of and to the west of the dense 
ridge; or the density of the molecular cloud itself 
must increase to the east, thereby inhibiting the 
expansion of the H II region on that side (Barnes 
et al. 1989). The contours of Figure 4 show that 
the dust extends further to the east of the ridge 
than to the west. If the source of the ionization 
front is to the west of the molecular ridge, then 
the vectors seen from our line of sight will stretch 
across the front of the ridge and to the east nearly 
in the plane of the sky (i.e., the Bios component 
is zero). To the west of the molecular ridge, such 
a source would produce an ionization front which 
would be bending around the ridge, with compo- 
nents Bios and B±_ both non-zero. 

However, there is substantial evidence that the 
dense cores and high density gas lie behind the H II 
region, as shown in Figure 6. Although Barnes et 
al. (1989) calculated the ionizing radiation de- 
tected in the nebula to require the IRS 3 source 
to be an 09 (or earlier) star, they conclude this 
to be unlikely unless the extinction toward IRS 3 
is 4.1 mag in a factor of 2 greater than the 
average extinction derived toward that star. In- 
stead, they suggest that an ensemble of B stars 
must have formed the H II region. If each star has 
formed a bubble of ionized radiation around it, as 
IRS 2 has produced the Eastern Loop (see Figure 
6), then the cumulative effect could produce the 



H II region. 

In order for this picture to be physically rea- 
sonable, the ionization front must be interacting 
with the molecular ridge. Gaume et al. (1992) 
detected 1.3 cm emission at 3" resolution located 
at south of FIR 4 (their NH3 4) and north of FIR 
5 (their NH3 5). North of the FIR 5 source, the 
1.3 cm emission arises from a dense, ionized in- 
terface between the molecular cloud (and ridge of 
cores) and the H II region. Ionized emission near 
FIR 4 has a compact component associated with 
the dense ridge. Additionally, higher NII3 temper- 
atures, larger line widths and smaller abundances 
suggest that the FIR 4 and FIR 5 sources are near- 
est the ionization front of the H II region. Schulz 
et al. (1991) show that ionized gas in fact fills the 
bay in the molecular bar surrounding IRS 3 (see 
their figure 7), and that the ridge is effectively sep- 
arated into two subregions by the ionization front 
of the H II region. We proceed on the basis that 
the ridge existed prior to the formation of the H II 
region. 

The polarization model shown in Figure 9 is 
based on the idea that the magnetic field has been 
swept up in a shell during the expansion of the 
H II region and stretched over the surface of the 
dense filamentary system of cores embedded in the 
wall of the molecular cloud. The stretching of the 
field in this scenario would tend to make the field 
quite uniform by the time the shell reached the 
system of cores. The structure of the magnetic 
field is idealized as follows. We assume that the 
field is swept into a sheet with a Gaussian planar 
distribution of the form 

p = p^e (7) 



when it first encounters the embedded filament, 
where y is along the line of sight, yc is the offset 
of the field's center from the filament axis, by is 
the width of the field region and B^^q and B^^^ 
are respectively the maximum, central values of 
the X and z field components. The magnetic field 
threading the shell is generally not aligned with 
respect to the filament axis, which is oriented in 
the z direction. The filament is assumed to be 
an unmagnetized, Ostriker (1964) filament, which 
we have truncated at a radius of 10*-^ tq, where rg 
is the radial scale of the filament given in Fiege & 
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Pudritz (2000a). This radial scale basically defines 
the radius where the inner, nearly flat density pro- 
file begins to steepen toward the profile that 
the density obeys at r ^ tq. 

We wrap the field threading the shell around 
the filament by applying a coordinate transforma- 
tion of the form 

x = xo, y = yo + ro{x/wf, z = zo (9) 

to the planar field giving by equation (8). The 

front of the shell is defined by the surface yo = 
— lO'-^ro, so that it just touches the filament. 
This transforms the originally planar sheet into 
a parabolic shell, wrapped around the filament, 
which is centered on the origin. The magnetic 
field and density are self-consistently distorted us- 
ing the Lagrangian forms of the MHD induction 
equation and the equation of mass conservation, 
as given by equation 6. The determinant |J| of 
the transformation is unity, so that the density is 
unchanged. 

The Zeeman data of Crutcher et al. (1999) 
suggest that the magnetic field is approximately 
parallel to the surface of the molecular cloud east 
of the filament, and strongly bent into the line of 
sight on the west. This is also supported by Barnes 
et al. (1989), who suggest that the molecular 
cloud may be denser east of the ridge than to the 
west, since the field would be swept further back 
on the side of lower density. Thus we rotate our 
field configuration by 20° counter-clockwise (look- 
ing toward the origin from the z axis) to mimic 
this asymmetry. The inclination angle that best 
reproduces the observed polarization map is 60° 
into the plane of the sky. We apply a final rota- 
tion of —15° in the plane of the sky to better agree 
with the observed orientation of the filament. The 
resultant polarization pattern is shown in Figure 
9. This figure also shows a map of the density- 
weighted average of the line of sight component of 
the magnetic field (arbitrary units) for comparison 
with the Zeeman data. 

This model is in better agreement with the Zee- 
man data for NGC 2024 than the heUcal field 
model given in § 4.3.1. However, our magnetized 
shell model has several shortcomings. The best- 
fitting polarization pattern that can be obtained 
from this type of model, with the Zeeman data 
regarded as a constraint, is of noticeably poorer 
agreement than the helical field model. The shell 



model reproduces the general features of the map, 
but the helical field model does a better job ac- 
counting for the detailed structure. The shell 
model also predicts a dense ridge of emission on 
the west side of the filament, where the parabolic 
shell plunges into the plane of the sky. It is note- 
worthy that the filament in this model is unmag- 
netized, and therefore does not contribute to the 
polarized flux. Hence, it is necessary to put a con- 
siderable amount of mass into the shell, so that 
the polarization percentage, defined as the polar- 
ized flux divided by the total intensity, is suffi- 
ciently high across the filament to account for the 
observations. The central density of the shell is 
11.4% of the central density of the filament, and 
the shell's thickness is 10'-^/3. A side-effect of this 
dense shell is that the emission away from the fila- 
ment is non-negligible. Both the polarization and 
the intensity in Figure 9 are clipped at 5% of their 
respective maximum values. Nevertheless, polar- 
ized emission above these levels is found in the 
background, in contrast to the helical field model, 
where this clipping procedure eliminates the fiux 
away from the filament. 

The depolarization observed in NGC 2024 is 
easily understood in this model. The filament it- 
self does not contribute to the polarized flirx, since 
it is unmagnetized. Instead, a uniformly polar- 
ized foreground sheet produces the polarization. 
Where the total intensity is high (i.e., along the 
filament), the polarization percentage declines, as 
we observe. 

A limitation of this model is that the filament 
is unmagnetized, and therefore obeys the den- 
sity profile predicted by Ostrikcr (1964). This den- 
sity gradient is known to be too steep to account 
for the observations in other filaments (Alves et 
al. (1998): Lada ct al. (1999): Johnstone & Bally 
(1999)). A magnetic field threading the filament 
could be included in this model, in principle, but 
this would make it substantially more complex. 

A third type of model may also be worth consid- 
ering. We have assumed in the above model that a 
shell swept up by the H II region wrapped around 
a dense filament, pre-existing in the cloud behind 
the H II region. One could imagine an alternate 
scenario involving the filament fragmenting out of 
the shell. The fastest growing unstable modes of 
fragmentation of a magnetized sheet result in file- 
ments oriented either parallel or perpendicular to 
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the magnetic field, depending on the thickness of 
the shell (Nagai et al. 1998). In either case, the 
geometry of the field probably could not account 
for the complex structure seen in our polarization 
map. We also note that a fragmenting shell would 
most likely produce more than one filament, and 
no others are observed in the vicinity of NGC 2024. 

4.4. Further Observational Tests of the 
Helical Model 

Although quantitative comparisons of the fila- 
mentary helically-threaded models to 850 /xm data 
(and where applicable, Zeeman data) will con- 
strain which model can reliably reproduce our ob- 
servations, different observation tests will even- 
tually be possible. For example, observations a) 
on larger spatial scales, b) over a range of wave- 
lengths, or c) at higher resolutions will each test 
the helical field models. 

If filaments arc threaded by helical fields, then 
the poloidal currents associated with those fields 
must connect to the larger scale magnetic field of 
the molecular cloud itself. The polarization sig- 
natures on these larger scales will be detectable 
once more sensitive submillimeter cameras (e.g. 
SCUBA-2) are built. Higher scnstitivity is only 
one goal of the next generation of detectors; the 
most important advance will come with the elim- 
ination of the need to chop on the sky. This 
will enable observations of polarized emission 
from faint, relatively featureless dust structures 
on scales larger than the detector field of view. 

The second extended test of helical field geome- 
tries can be made with high resolution, interfero- 
metric, mapping. The helical field geometry pre- 
dicts that depolarization should persist down to 
very small scales. In practice, such maps are hard 
to interpret since current interferometers are lim- 
ited to observations of cores where star formation 
is already occurring. Ideally, one would like to ob- 
serve filaments at equally high resolution, but this 
goal is hindered by sensitivity limits and the fil- 
tering of large scale structure beyond some limit 
imposed by the spacing of the array. While sensi- 
tivity will continue to improve, the filtering effect 
is not so easily overcome where short spacings in- 
formation is not present. 

Perhaps the most immediately tangible way 
to further test helical field models is the obser- 



vation of regions over a range of wavelengths. 
If one were to observe shorter wavelengths until 
sources become optically thick, then the measure- 
ments could be used as a depth probe of the mag- 
netic field. Although grain effects can produce 
changes in polarization percentage (Hildebrand et 
al. 1999), changes in position angle are attributed 
to the magnetic field. If a field is indeed heli- 
cal, the vectors should change orientation relative 
to a straight filament, from alignment with the 
filament axis to orthogonality as wavelength de- 
creases. The future SOFIA polarimeter, operat- 
ing below 200 /im will be particularly useful for 
testing this aspect of the models. 

5. Summary 

The polarization patterns observed in Orion B 
are dominated by orderly structure on scales at 
least as large as the areas mapped (with the possi- 
ble exception of the LBS 23N region). This rcsidt 
can only be consistent with magnetic fields that 
are ordered on similar scales or larger. The fila- 
mentary clouds, NGC 2024 and LBS 23N, exhibit 
polarization patterns that are inconsistent with 
purely poloidal magnetic fields threading the fila- 
ments, or with fields simply threaded transversely 
to them. Curved field lines are necessary to model 
the polarization data in NGC 2024. The fact that 
the polarization pattern is symmetric about the 
dense ridge of cores suggests that there is a corre- 
lation between the presence of dense gas and the 
exhibited polarization pattern. 

The polarization systematically decreases with 
intensity for all three regions mapped. This result 
is in agreement with Paper II, which notes depo- 
larization toward the axis of the Integral-shaped 
Filament in OMC-3, as well as most submillime- 
ter and millimeter observations of star-forming re- 
gions (see Weintraub, Goodman & Akeson (2000) 
for a review). This "polarization hole" effect is a 
global property of our maps, and not just a local 
phenomenon that occurs for a few pixels near the 
very brightest emission. 

NGC 2071 is a massive core forming multiple 

protostars. Its polarization pattern is ordered and 
qualitatively similar to that of OMC-1, which is 
an even more massive core in Orion A. In OMC- 
1, Schlciming (1998) interpreted their polarization 
data as being due to an hourglass magnetic field, 
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with the field lines pinched due to the collapse of 
the core. This interpretation is consistent with the 
flattening observed along the inferred field lines in 
OMC-1. On the other hand, NGC 2071 docs not 
show any flattening, which is inconsistent with an 
interpretation of the magnetic field threading this 
core as a dynamically significant hourglass field. 
We note that if the vectors of NGC 2071 were ro- 
tated to infer a net field direction, this direction 
would be perpendicular to the most powerful out- 
flow in the region. Alternatively, our map could 
be interpreted as resulting from a field that is pre- 
dominantly toroidal about the axis of symmetry. 
Although CO J = 3 — 2 emission can contam- 
inate SCUBA 850 nm data, we estimate it can 
contribute a maximum of 1% of polarized emission 
when optically thick. Since much higher polariza- 
tions are measured, we conclude that the polar- 
izations measured in NGC 2071 arc dominated by 
dust emission. Application of the Chandrasekhar 
& Fermi method to this core yields a mean field 
strength estimate of 56 fiG, and we find that half 
the magnetic energy is accounted for by the mean 
field component. 

LBS 23N exhibits the least orderly polarization 
patterns of the three regions studied. Most of the 
polarization vectors are aligned in a north-south 
orientation, particularly to the east of the cores. 
The cores themselves are significantly depolarized. 
Near the southern boundary of the map, the vec- 
tors rotate by 90° to an east-west orientation. This 
abrupt change might be explained either by an ex- 
tension of the Ficgc & Pudritz model (sec §4.2), or 
by a smaller filament orthogonal to the main fila- 
ment. More data to the south of LBS 23N would 
be needed to distinguish between these possibili- 
ties. 

The NGC 2024 ridge of cores presents the most 
interesting polarization data of this paper. The 
polarization is strong everywhere except the cores, 
which are significantly depolarized. Our map 
agrees very well with the 100 /im polarization data 
obtained with the KAO Stokes polarimeter (Dot- 
son et al. 2000), even though these wavelengths 
likely probe different dust temperatures. Interest- 
ingly, the region between FIR 4 and FIR 5 ex- 
hibits diminished intensity but no polarization is 
detected there. We are able to successfully model 
the polarization pattern using a helical field geom- 
etry threading a curved filament. However, given 



the overall geometry of the region an H II region 
expanding into the molecular cloud in which the 
cores are embedded - we slightly prefer a model in 
which the ionization front has swept up the mag- 
netic field and is now stretching it around the ridge 
of dense cores. Since this model contains an un- 
polarized filament, it may be oversimplified and a 
more complete, complex picture will include both 
the polarized ionization front and a polarized fila- 
ment. 

We note that our analysis of the NGC 2024 po- 
larization data have been aided immensely by the 
existence of Zeeman maps of the line-of-sight field 
strength toward part of this region published by 
Crutcher et al. (1999). Unfortunately, very few 
regions have been successfully mapped using Zee- 
man splitting (Crutcher 1999). This is an area 
where future technological advances in instrumen- 
tation would greatly improve our understanding of 
magnetic fields in molecular clouds. 
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A. Polarization Data 



Tables 1, 2, and 3 present the polarization data for the regions NGC 2071, LBS 23N, and NGC 2024 
respectively. Each table contains only those vectors plotted on Figures 1, 2, and 3. Positional offsets 
for NGC 2071 are from aj20oo = 05'^47™5n6 (05'^44™31!0 in B1950) and ^j2ooo = +00° 21' 47'/l (+00° 
20' 45'.'0 in B1950). The reference position for LBS 23N is ajaooo 05h46'"ll!16 (05'^43'"37r0 in B1950) 
and 5j2000 = -00° 09' IS'.'O (-00° 10' 15'.'2 in B1950). Finally the reference position for NGC 2024 is 
Q!j20oo = 05Ml'"44!l (05'^39™12!6 in B1950) and 5j20oo = -01° 54' 49'.'7 (-01° 56' 15'.'0 in B1950). 

The columns in each table are: (1) the R.A. offset in arcseconds; (2) the DEC. offset in arcseconds; (3) the 
polarization percentage at that position; (4) the uncertainty in the polarization percentage; (5) the signal- 
to-noise ratio of the polarization percentage; (6) the polarization position angle; and (7) the uncertainty in 
the polarization position angle. 
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Table 1 

NGC 2071 850 //M Polarization Data 



A R.A. A DEC. p 

n (") (%) 



-16.5 


-67.5 


5.56 


10.5 


-58.5 


12.99 


1.5 


-58.5 


8.51 


-7.5 


-58.5 


5.14 


-25.5 


-58.5 


6.04 


-43.5 


-58.5 


6.64 


10.5 


-49.5 


10.04 


1.5 


-49.5 


5.25 


-25.5 


-49.5 


4.73 


-43.5 


-49.5 


13.76 


-52.5 


-49.5 


7.48 


10.5 


-40.5 


3.54 


-16.5 


-40.5 


8.00 


-25.5 


-40.5 


4.56 


-34.5 


-40.5 


9.33 


-52.5 


-40.5 


15.03 


37.5 


-31.5 


9.70 


28.5 


-31.5 


11.90 


19.5 


-31.5 


3.11 


10.5 


-31.5 


2.49 


-16.5 


-31.5 


2.73 


46.5 


-22.5 


4.34 


19.5 


-22.5 


1.77 


10.5 


-22.5 


2.48 


1.5 


-22.5 


1.25 


-7.5 


-22.5 


1.83 


-16.5 


-22.5 


1.66 


-34.5 


-22.5 


5.31 


-52.5 


-22.5 


8.94 


37.5 


-13.5 


5.70 


28.5 


-13.5 


2.46 


19.5 


-13.5 


2.29 


10.5 


-13.5 


1.28 


1.5 


-13.5 


1.04 


-7.5 


-13.5 


1.53 


-34.5 


-13.5 


3.15 


-43.5 


-13.5 


5.76 


-52.5 


-13.5 


4.17 


37.5 


-4.5 


4.77 


10.5 


-4.5 


1.74 


1.5 


-4.5 


1.09 



d/p 


Up 


\ ) 


d6 
\ ) 


1.07 


5.2 


63.8 


5.5 


1.41 


9.2 


-14.7 


3.1 


1.33 


6.4 


7.0 


4.5 


1.04 


4.9 


42.8 


5.8 


1.40 


4.3 


-47.7 


6.7 


1.31 


5.1 


-20.4 


5.7 


1.31 


7.7 


-29.3 


3.7 


0.90 


5.8 


18.8 


4.9 


0.87 


5.4 


-80.1 


5.3 


1.48 


9.3 


-8.0 


3.1 


1.29 


5.8 


22.0 


5.0 


0.61 


5.9 


-16.1 


4.9 


0.48 


16.6 


-14.2 


1.7 


0.68 


6.7 


6.9 


4.3 


1.46 


6.4 


-19.2 


4.5 


1.40 


10.7 


-17.2 


2.7 


1.41 


6.9 


46.9 


4.2 


1.02 


11.6 


26.4 


2.5 


0.59 


5.2 


-2.2 


5.5 


0.35 


7.1 


13.0 


4.1 


0.35 


7.9 


47.0 


3.6 


1.06 


4.1 


54.2 


7.0 


0.34 


5.2 


-1.2 


5.5 


0.23 


10.7 


26.1 


2.7 


0.26 


4.8 


47.0 


6.0 


0.21 


8.8 


70.0 


3.2 


0.25 


6.6 


-85.7 


4.4 


0.88 


6.0 


21.6 


4.7 


0.83 


10.8 


49.6 


2.7 


0.93 


6.1 


57.6 


4.7 


0.57 


4.3 


37.1 


6.6 


0.30 


7.6 


13.5 


3.8 


0.21 


6.0 


18.6 


4.7 


0.17 


6.3 


35.0 


4.5 


0.18 


8.4 


53.0 


3.4 


0.31 


10.1 


57.8 


2.8 


0.50 


11.6 


37.1 


2.5 


0.68 


6.1 


30.2 


4.7 


0.97 


4.9 


18.0 


5.8 


0.17 


10.3 


14.1 


2.8 


0.10 


11.2 


29.7 


2.5 
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Table 1 — Continued 



A R.A. A DEC. p 

n (") (%) 



-7.5 


-4.5 


1.14 


-25.5 


-4.5 


1.20 


-34.5 


-4.5 


2.64 


-43.5 


-4.5 


6.45 


-52.5 


-4.5 


7.21 


37.5 


4.5 


8.08 


28.5 


4.5 


2.78 


-7.5 


4.5 


1.24 


-16.5 


4.5 


1.66 


-25.5 


4.5 


1.44 


-43.5 


4.5 


3.73 


37.5 


13.5 


5.95 


28.5 


13.5 


1.43 


19.5 


13.5 


3.03 


10.5 


13.5 


1.31 


-7.5 


13.5 


1.41 


-16.5 


13.5 


1.47 


-25.5 


13.5 


2.18 


-34.5 


13.5 


5.23 


-43.5 


13.5 


4.95 


-52.5 


13.5 


5.38 


46.5 


22.5 


8.81 


37.5 


22.5 


3.78 


19.5 


22.5 


2.73 


10.5 


22.5 


1.31 


1.5 


22.5 


2.42 


-7.5 


22.5 


1.64 


-16.5 


22.5 


2.35 


-25.5 


22.5 


2.60 


-34.5 


22.5 


3.45 


-43.5 


22.5 


5.42 


-52.5 


22.5 


9.78 


46.5 


31.5 


11.95 


37.5 


31.5 


6.91 


19.5 


31.5 


4.49 


10.5 


31.5 


4.60 


-25.5 


31.5 


2.74 


-34.5 


31.5 


4.45 


-43.5 


31.5 


6.17 


37.5 


40.5 


5.39 


19.5 


40.5 


4.60 



dp 


dp 





d0 


(%) 




(°) 


(°) 


0.09 


13.3 


55.8 


2.2 


0.19 


6.2 


55.5 


4.6 


0.26 


10.2 


75.3 


2.8 


0.60 


10.7 


37.7 


2.7 


0.69 


10.4 


13.9 


2.8 


0.81 


10.0 


24.3 


2.9 


0.39 


7.2 


-20.5 


4.0 


0.15 


8.3 


61.8 


3.4 


0.17 


9.5 


54.2 


3.0 


0.22 


6.7 


44.0 


4.3 


0.59 


6.3 


37.1 


4.5 


0.76 


7.8 


21.6 


3.7 


0.33 


4.3 


-15.0 


6.7 


0.25 


12.4 


-2.0 


2.3 


0.25 


5.3 


-2.7 


5.4 


0.24 


5.9 


50.4 


4.9 


0.20 


7.3 


58.9 


3.9 


0.32 


6.8 


38.0 


4.2 


0.44 


12.0 


21.3 


2.4 


0.57 


8.6 


19.3 


3.3 


0.84 


6.4 


16.8 


4.5 


1.23 


7.2 


28.4 


4.0 


0.65 


5.8 


25.3 


4.9 


0.29 


9.3 


29.7 


3.1 


0.29 


4.5 


26.9 


6.4 


0.32 


7.4 


74.1 


3.8 


0.32 


5.0 


60.2 


5.7 


0.33 


7.1 


61.0 


4.0 


0.34 


7.6 


31.4 


3.8 


0.55 


6.3 


21.2 


4.5 


0.54 


10.0 


17.7 


2.9 


0.68 


14.4 


17.8 


2.0 


1.13 


10.6 


27.9 


2.7 


0.72 


9.6 


26.4 


3.0 


0.86 


5.2 


52.1 


5.5 


0.80 


5.8 


17.6 


5.0 


0.64 


4.3 


-0.8 


6.7 


0.51 


8.8 


13.6 


3.3 


0.83 


7.5 


-7.8 


3.8 


0.73 


7.4 


-84.2 


3.9 


0.83 


5.5 


0.8 


5.2 
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Table 1 — Continued 



A R.A. 

(") 


A DEC. 

(") 


P 

(%) 


dp 
(%) 


CTj, 


e 

n 


de 

n 


1.5 


40.5 


6.54 


0.79 


8.2 


14.3 


3.5 


-34.5 


40.5 


3.77 


0.78 


4.8 


43.7 


5.9 


-43.5 


40.5 


9.15 


0.97 


9.4 


17.7 


3.0 


1.5 


49.5 


8.38 


1.12 


7.5 


17.0 


3.8 


-7.5 


49.5 


7.30 


1.04 


7.0 


14.0 


4.1 


-16.5 


49.5 


9.14 


0.96 


9.5 


50.5 


3.0 


-34.5 


49.5 


6.52 


0.84 


7.8 


80.0 


3.7 


37.5 


58.5 


3.30 


0.80 


4.1 


-38.6 


7.0 


28.5 


58.5 


5.63 


0.71 


7.9 


15.2 


3.6 


19.5 


58.5 


4.63 


0.73 


6.3 


25.0 


4.5 


1.5 


58.5 


6.19 


1.22 


5.1 


18.3 


5.7 


-7.5 


58.5 


5.60 


1.31 


4.3 


8.9 


6.7 


-16.5 


58.5 


8.65 


1.01 


8.5 


19.5 


3.4 


-25.5 


58.5 


3.16 


0.69 


4.5 


35.9 


6.3 


-34.5 


58.5 


3.84 


0.91 


4.2 


33.1 


6.8 


19.5 


67.5 


14.62 


1.33 


11.0 


64.9 


2.6 


10.5 


67.5 


10.32 


1.19 


8.7 


36.5 


3.3 
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Table 2 

LBS 2:-)N 850 /m Polarization Data 



A R.A. A DEC. p 

n (") (%) 



-9.0 


-124.5 


5.58 


15.0 


-112.5 


12.03 


3.0 


-112.5 


4.29 


-9.0 


-112.5 


1.83 


15.0 


-100.5 


11.75 


3.0 


-100.5 


5.34 


-9.0 


-100.5 


3.14 


-21.0 


-100.5 


2.88 


-33.0 


-100.5 


3.12 


-45.0 


-100.5 


8.47 


15.0 


-88.5 


8.27 


3.0 


-88.5 


4.88 


-9.0 


-88.5 


1.97 


-21.0 


-88.5 


2.94 


-33.0 


-88.5 


3.77 


-45.0 


-88.5 


5.28 


27.0 


-76.5 


13.00 


15.0 


-76.5 


5.45 


3.0 


-76.5 


3.61 


-9.0 


-76.5 


6.55 


-21.0 


-76.5 


7.21 


-45.0 


-76.5 


12.16 


27.0 


-64.5 


7.82 


15.0 


-64.5 


3.21 


-9.0 


-64.5 


2.52 


-21.0 


-64.5 


3.40 


27.0 


-52.5 


12.61 


15.0 


-52.5 


6.77 


3.0 


-52.5 


1.47 


-9.0 


-52.5 


2.61 


-21.0 


-52.5 


2.09 


-33.0 


-52.5 


3.51 


27.0 


-40.5 


8.82 


15.0 


-40.5 


2.61 


3.0 


-40.5 


1.55 


27.0 


-28.5 


6.18 


15.0 


-28.5 


5.57 


3.0 


-28.5 


1.58 


-9.0 


-28.5 


1.40 


-21.0 


-28.5 


1.98 


-33.0 


-28.5 


5.60 



d/p 


Up 


\ ) 


d6 
\ ) 


1.03 


5.4 


-29.1 


5.3 


1.37 


8.8 


29.9 


3.3 


0.60 


7.1 


38.4 


4.0 


0.30 


6.1 


-16.7 


4.7 


1.25 


9.4 


5.1 


3.1 


0.68 


7.8 


74.4 


3.7 


0.34 


9.1 


-82.0 


3.1 


0.30 


9.7 


-88.8 


3.0 


0.62 


5.0 


71.7 


5.7 


1.35 


6.3 


75.1 


4.5 


1.12 


7.4 


-11.9 


3.9 


0.68 


7.1 


-3.4 


4.0 


0.49 


4.1 


-88.2 


7.1 


0.50 


5.9 


29.4 


4.9 


0.76 


5.0 


71.7 


5.7 


1.28 


4.1 


51.7 


6.9 


1.34 


9.7 


28.8 


3.0 


0.80 


6.9 


21.6 


4.2 


0.65 


5.6 


-26.6 


5.1 


0.57 


11.5 


-53.6 


2.5 


0.77 


9.4 


-42.3 


3.1 


1.20 


10.2 


7.8 


2.8 


1.29 


6.0 


1.8 


4.7 


0.63 


5.1 


10.3 


5.6 


0.47 


5.4 


-78.8 


5.4 


0.71 


4.8 


62.8 


5.9 


1.01 


12.5 


8.1 


2.3 


0.47 


14.5 


4.9 


2.0 


0.26 


5.7 


16.3 


5.0 


0.23 


11.5 


5.2 


2.5 


0.42 


5.0 


50.1 


5.7 


0.67 


5.3 


0.6 


5.4 


0.74 


11.9 


6.9 


2.4 


0.37 


7.1 


32.3 


4.1 


0.13 


11.6 


23.5 


2.5 


1.06 


5.8 


38.5 


4.9 


0.41 


13.5 


28.3 


2.1 


0.18 


8.8 


27.8 


3.2 


0.13 


10.9 


46.0 


2.6 


0.24 


8.4 


7.9 


3.4 


0.63 


8.9 


-6.9 


3.2 
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Table 2 — Continued 



ADA 

A H.A. 




P 


dp 






dv 


( ) 


( ) 


(%) 


(%) 


( ) 


( ) 


51.0 


-16.5 


8.44 


1.14 


7.4 


17.0 


3.9 


15.0 


-16.5 


7.53 


0.92 


8.2 


29.0 


3.5 


3.0 


-16.5 


2.23 


0.41 


5.5 


11.5 


5.2 


-9.0 


-16.5 


1.47 


0.30 


4.9 


36.3 


5.9 


-21.0 


-16.5 


2.83 


0.36 


7.9 


15.6 


3.6 


-45.0 


-16.5 


12.59 


1.23 


10.2 


30.8 


2.8 


39.0 


-4.5 


6.46 


1.17 


5.5 


-4.9 


5.2 


15.0 


-4.5 


21.96 


1.30 


16.8 


-4.4 


1.7 


3.0 


-4.5 


2.54 


0.54 


4.7 


15.1 


6.1 


-33.0 


-4.5 


4.97 


0.58 


8.6 


-29.4 


3.3 


39.0 


7.5 


4.09 


0.77 


5.3 


29.8 


5.4 


27.0 


7.5 


8.46 


1.11 


7.7 


19.7 


3.7 


15.0 


7.5 


6.71 


0.83 


8.1 


16.4 


3.6 


3.0 


7.5 


4.21 


0.53 


8.0 


-28.8 


3.6 


-9.0 


7.5 


1.17 


0.23 


5.0 


45.7 


5.8 


-21.0 


7.5 


1.26 


0.28 


4.5 


-10.5 


6.3 


-33.0 


7.5 


3.28 


0.50 


6.5 


20.0 


4.4 


39.0 


19.5 


11.24 


1.09 


10.3 


-62.0 


2.8 


27.0 


19.5 


10.48 


1.05 


9.9 


-30.8 


2.9 


15.0 


19.5 


5.70 


0.58 


9.8 


21.4 


2.9 


3.0 


19.5 


4.01 


0.44 


9.2 


-23.6 


3.1 


-9.0 


19.5 


5.57 


0.43 


12.9 


14.6 


2.2 


-21.0 


19.5 


2.68 


0.44 


6.1 


14.9 


4.7 


-33.0 


19.5 


7.10 


0.43 


16.4 


-6.7 


1.8 


-45.0 


19.5 


1.89 


0.46 


4.1 


26.3 


7.0 


3.0 


31.5 


5.46 


0.66 


8.2 


57.2 


3.5 


-21.0 


31.5 


11.41 


0.66 


17.3 


-4.8 


1.7 


-33.0 


31.5 


5.90 


0.57 


10.3 


36.8 


2.8 


-45.0 


31.5 


7.10 


0.78 


9.1 


11.6 


3.2 
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Table 3 

NGC 2024 850 //M Polarization Data 



A R.A. 


A DEC. 


P 


dp 


CTp 


e 


de 


(") 


(") 


(%) 


(%) 


n 


n 
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Fig. 1.— 850 fim polarization of NGC 2071. NGC 2071IR greyscale from (Mitchell et al. 2001) overlain 
with polarization vectors. The vectors are binned to 3 pixels (9") so each vector is separated by just over 
half the 850 //m beamwidth (14" at the JCMT). Vectors are only plotted where uncalibrated / > 2% of 
the source peak, polarization percentage exceeds 1% and the absolute uncertainty in polarization percentage 
is no greater than 1.5%. All vectors plotted have signal-to-noise in polarization percentage, Cp > 4 (and 
therefore an uncertainty in position angle, dO < 7.2°), but bold vectors have CTp > 6 {d6 < 4.8°). Over the 
entire region, the mean position angle is 22° while the standard deviation is 32°. The mean polarization 
percentage is 5.1% with a standard deviation of 3.3%. Cores are labelled as in Mitchell et al. (2001), and 
blue lines mark the mean direction of the outflow from IRS 3. Right: Details of the central core (# 11 of 
Mitchell et al. (2001)) with the same vector cutoffs except / > 12% of the peak flux. Blue contours trace 
the surface density in unpolarized 850 fim intensity. All coordinates are J2000. 
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Fig. 2.— 850 /xm polarization of LBS 23N. LBS 23 greyscale is taken from Mitchell et al. (2001). The 
polarization vectors overlain have been binned to 12". Vectors are plotted where / > 4% of the peak of core 
59, p > 1%, dp < 1.5% and p/dp > 4 {d0 < 7.2°). Bold vectors show those where p/dp > 6 {d0 < 4.8°). 
The core identifications are taken from Mitchell et al. (2001). The mean polarization percentage measured 
is 5.6% with a standard deviation of 3.8%. A potential "cross" filament could be oriented east-west and has 
been indicated by the presence of a dashed line. Coordinates are J2000. 
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Fig. 3. — NGC 2024 850 polarization. NGC 2024 greyscale overlain with polarization 12" binned 
polarization vectors, almost the JCMT beamwidth of 14" at 850 /im. Vectors are plotted where / > 2.5% 
of the FIR 5 peak, p > 1%, dp < 1% and p/dp > 6 {d6 < 4.8°). Bold vectors show those where p/dp > 10 
{dO < 2.9°). Overall, these data represent the highest signal to noise polarization detections in our data 
set. The mean percentage polarization is 3.4% with a standard deviation of 2.3% in 159 vectors. The 
submillimeter cores are labelled as in Mezger et al. (1988). We have indicated the orientation of the FIR 
5 unipolar outflow by a blue dashed line, and the CO J = 3 — 2 line polarization vectors have been plotted 
(in blue) as measured by Greaves et al. (2001a). Coordinates are J2000. 
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Fig. 4. — Comparison of Polarizations at 100 /xm and 850 /im. Comparison of the polarization patterns at 
850 fim (red) and 100 fim (green) reveals a strong consistency in position angle. The polarization percentages 
are considerably lower for 100 /im, possibly indicating that fewer dust grains are sampled within hotter dust 
or, more probably, that the dust grain alignment is lessened at higher temperatures. The 850 /^m contours 
are plotted at 80, 90, 95, 98, 99, 99.5 and 99.9 percentiles. A blue cross marks the location of the maximum 
line-of-sight magnetic field as measured by Crutcher et al. (1999). 
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Fig. 5. - The Depolarization Effect in Orion B. The depolarization effect typically seen in protostars and 
extended regions (such as OMC-3 in Orion A) is also observed in NGC 2071, NGC 2024 and LBS 23N. 
For all three regions, the brightest positions have the smallest polarizations, exhibiting < 2% whereas the 
faintest regions show polarizations approaching and exceeding 10%. Decreased polarization percentage is 
observed for a large range of / values in NGC 2024. This region's cores also contain the largest areas where 
depolarization is total, or the levels of polarization are so low (i.e., < 1%) as to be consistent with zero. 
Extremely high volume densities (on the order of 10^ cm~^) have been estimated in NGC 2024 (Mezger et 
al. 1988). 
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Fig. 6. — Schematic Geometry of the NGC 2024 region. The geometry of the NGC 2024 H II dense cores 
and foreground dark dust lane, as seen from the western edge of the cloud. From the' Earth's perspective, 
the dark lane (unassociated with the dense cores) is seen in the foreground of the H II region. The H II 
region has blistered out of the molecular cloud to the north, west and east but is hindered to the south by 
dense molecular material as indicated by radio maps of the ionized emission (Barnes et al. 1989). The IR 
sources marked include IRS 2 (within the Eastern Loop) and IRS 3, near the southern boundary of the H II 
region. 
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Fig. 7. — A Helical Field Model for NGC 2024. We show an example of a "bent filament" model (discussed 
in Section 4.3) for comparison with NGC 2024. The model is a Fiege & Pudritz (2000a) filament model 
constructed with the following parameters: C = 1.1, Vz = 8, F^ = 15. The length of the filament is 6 x 10*-^ 
and the ends have been rounded. We have bent the entire filament into a circular arc perpendicular to the 
plane of the sky and toward the observer, keeping the top of the filament parallel to the original orientation. 
The same pattern would be produced for a filament bent away from the observer. The radius of the arc is 
3/7r times the filament length. We then rotated the entire structure by 20° and inclined it relative to the 
plane of the sky by —15°. The right panel shows the expected line-of-sight field pattern from a helical field. 
A reversal in field direction is expected across the filament. 
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Fig. 8. — Geometry of the NGC 2024 region, viewed from the north. The geometry of the NGC 2024 region 
as viewed from the north, looking down on the ridge of dense cores. A proposed orientation of the field line 
is shown by the dot-dashed line. We speculate that the field has been compressed into a uniform layer on 
the edge of the expanding H II region. The ionization front has penetrated the molecular material more 
deeply to the west where the cloud is less dense, but is hindered by the dense ridge. This accounts for the 
field lying predominantly in the plane of the sky to the east of the ridge, but significantly in the line of sight 
to the west. 
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Fig. 9. — Polarization Pattern Produced by a Parabolic Field. Model showing the polarization pattern from 
a parabolic field bent around a dense ridge of gas, under the scenario depicted in Figure 8 and described 
in the text. The left panel shows the polarization pattern and the right traces the Bios- The filament is 
unmagnetized and the field is contained entirely in the shell of gas. The density of the shell containing the 
field must be ^ 10% that of the central axis of the filament. The parabolic shell is originally symmetric 
along the line of sight, but is rotated by 20° to produce the correct behavior in the polarization and Bios- 
The inclination is 60° and the whole filament is rotated by on the plane of the sky —15°. 
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